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Life after eruption - II. The eclipsing old nova V728 Scorpii 
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ABSTRACT 

The old nova V728 Sco has been recently recovered via photometric and spectroscopic obser- 
vations, 150 years after the nova eruption. The spectral properties pointed to a high-inclination 
system with a comparatively low mass-transfer rate. In this paper we show that the object is 
an eclipsing system with an orbital period of 3.32 h. It has enhanced long-term variability 
that can be interpreted as 'stunted' dwarf-nova-type outbursts. Using the ingress and egress 
times of the eclipsed components we calculate the radius of the central object. The latter turns 
out to be significantly larger than a white dwarf and we identify it with a hot inner disc. The 
implications for models on the behaviour of post-novae are discussed. 

Key words: binaries: close - eclipsing - novae, cataclysmic variables 



1 INTRODUCTION 

A nova eruption is a thermonuclear explosion on the surface of 
an accreting white dwarf. In the process, a part of the outer lay- 
ers of the white dwarf is ejected as a nova shell. It is not clear if 
the expelled material amounts to less or more than the accreted 
mat erial, although the newest studies point t o the former possibil- 
ity dZorotovic. Schreiber. & Gansickell20l"li) . In the case of clas- 
sical novae, the donor usually is a late-type main- sequence star. 
This makes the system hosting the nova event a cataclysmic vari- 
able (CV). Since the CV is not destroyed by the explosion, it is 
safe to assume that a nova eruption is a recurrent even t. A typical 
shell mass of 10~^-10"^ M© fe.g./Yaron eTaDl2005h and an up- 
per limit of the average mass-transfer rate of ~ 10~^ Mq yr~^ 
jTownslev & Gansickd l2009t) then translates to recur rence times 
> 10^ yr. This is now also supported observationally (IShara et al.l 
l2012ah . 

Nova eruptions can thus be expected to be part of the evolution 
of every CV, assuming that the mass-transfer rate is high enough to 
accumulate the necessary critical amount of material on the surface 
of the white dwarf within the lifetime of the system. Hence, once 
the effects of the nova eruption have worn off, the behaviour of the 
underlying CV should be determined (once again) by the properties 
that make it part of one of the many subtypes of CVs, i.e. especially 
its orbital period, mass-transfer rate and the magnetic field strength 
of the white dwarf. There are several post-novae that corroborate 
this assumption, e.g. DQ Her (N ova Her 1934) is the prototype 
intermediate polar (iWa^lT983t) . and RR Pic (Nova Pic 1925) 



Table 1. Log of observations. 



Date 


Filter/grism (A) 


n 


^exp (s) 


At (h) 


2012-03-25 


#20 [6035-7135, 3.8] 


3 


900 


1.78 


2012-03-26 


#20 [6035-7135, 3.8] 


24 


600/300 


2.37 


2012-03-27 


#20 [6035-7135, 3.8] 


41 


300 


3.84 


2012-03-28 


V 


299 


20 


4.59 


2012-04-01 


V 


151 


20 


2.54 


2012-05-15 


V 


61 


40 


1.29 


2012-05-16 


#7 [3015-5190, 7.0] 


12 


300 


1.04 




#16 [5965-10250, 13.6] 


10 


300 


1.09 



belon gs to the SW Sex stars JSchmidtobreick. Tappert. & Saviand 
l2003h . Of special importance in this context is the discovery of 
an ancient nova shell around the dwarf novae Z Cam and AT Cnc 
( Shara et al .120071 , 120i2b), which made these the first CVs to be es- 
tablished as post-novae not based on the observation of the actual 
nova eruption. 

Within our project to establish a sample of post-novae to com- 
pare their characteristics as a group we have recently recovered the 
position of the system V728 Sco (Nova S co 1862; TebbuttlllSTSh 
150 yr after its eruption (iTappert et al. 12012. hereafter Paper I). The 
spectral characteristics indicated a high system inclination and a 
comparatively low mass-transfer rate. Here we present a more de- 
tailed study of that object. 



* E-mail: claus. tappert @uv.cl 

t Based on observations with ESO telescopes, proposal numbers 088.D- 
0588 and 089.D-0505 



2 OBSERVATIONS AND REDUCTION 



We observed V728 Sc o during observing runs in 2012 IMa rch and in 
May using EF0SC2 jEckert. Hofstadt. & ]V[eln"ic3l 19891) mounted 
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on the ESO-NTT at La Silla, Chile. We employed the Bessell V fil- 
ter to obtain a photometric light curve, the medium resolution grism 
#20 to measure radial velocities and the two low resolution grisms 
#7 and #16 to examine the spectral variation of the blue and the red 
part of the spectrum throughout eclipse. The observations are sum- 
marized in Table[T] It states the start of the night of the observations 
and the filter or grism used for the photometric or spectroscopic ob- 
servations, respectively. The column for the latter also includes the 
nominal wavelength range and the spectral resolution as the typi- 
cal full width at half-maximum (FWHM) of an arc line. The actual 
useful wavelength range is somewhat shorter for the low resolution 
grisms: due to the efficiency limits the signal-to-noise ratio (S/N) 
decreases strongly at wavelengths <3650 A (grism #7) and >9500 
A (grism #16). The latter grism additionally suffers from fringing 
for wavelengths >7000 A. Table [T] furthermore contains informa- 
tion on the number n of frames taken, the exposure time texp of 
an individual data frame, and the time range At covered by the 
observations. 



2.1 Photometry 

The reduction of the photometric data includes the subtraction of 
bias frames, but no flat-field correction because EF0SC2 flats are 
affected by a central light concentration. Using iRAF's DAOPHOT 
package, relative photometric magnitudes were obtained for all 
stars that were found by the DAOFIND algorithm. On some frames 
during deep eclipse, the target (V728 Sco) was not detected by that 
routine, and its position was determined with respect to a number 
of neighbour stars. An aperture radius of 3.0 pixels was used for 
all measurements. The data subsequently served a s input for th e 
standalone DAOMATCH and DAOMASTER routines dStetsonll 1992b . 
Differential magnitudes for V728 Sco were computed with respect 
to the average of 1 1 suitable comparison stars in the neighbourhood 
(dz200 pixels). Since no standard star data were taken during these 
observing runs, calibrated magnitudes w ere obtai ned by compari- 
son with the photometric data reported in Paper I. We furthermore 
computed differential light curves for all other stars on the field, but 
did not find any clear intrinsic variability for any of them. 
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2.2 Spectroscopy 

The spectroscopic data were reduced by subtracting the nightly av- 
erage bias frame and by dividing them by an average flat-field. 
The latter had been normalised for its spectral energy distribution 
by division through a high-order cubic spline fit along the disper- 
sion axis. The spectra were subsequently extracted using the im- 
plementation of the optimal extraction algorithm (Home 1986) in 
IRAF's APALL routine. Wavelength calibration was performed us- 
ing Thorium-Argon spectra that were taken during the afternoon. 
Flux-calibrated spectra for grisms #7 and #16 were obtained in 
comparison to the standard star LTT 7987 that was observed during 
the same night. Since an absolut flux calibration was not necessary 
for our purposes we did not perform any correction for slit losses. 
No flux standards were observed for the grism #20 spectra, which 
nevertheless were corrected for the instrumental response function 
using previously taken data. 



Figure 1. Orbital light curves from the photometric V data and from nor- 
malised fluxes of the spectroscopic data taken with grisms #20, #7, and #16. 
The data are folded on the ephemeris given in equation ([T] 



3 RESULTS 

3.1 Light curves and the orbital period 

The photometric data yield the light curves shown in the top three 
panels in Fig.[Tl The most prominent feature is doubtlessly the deep 
eclipse. We furthermore notice the presence of an orbital hump that 
lets the photometric brightness rise to a maximum just prior to the 
eclipse. This hump, which most probably corresponds to the hot 
spot, i.e. the region where the accretion from the secondary star 
impacts on the accretion disc, is more prominent in the second of 
the light curves, from April 1^^, than in the other two, from March 
28 and from May 15. In fact, the April data show the system in 
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Table 2. Long-term photometric data. 



HJD-2450000 


R (mag) 


V (mag) 


Notes 


4972.5 




18.47(11) 


Calibrated data 


5742.5 


17.431(03) 


18.03(01) 


Acquisition frame 


6012.5 


16.41(02) 


17.01(02) 


Acquisition frame 


6013.5 


16.49(03) 


17.08(03) 


Acquisition frame 


6014.5 


16.61(02) 


17.21(02) 


Acquisition frame 


6015.5 




17.33(09) 


Light curve 


6019.5 




18.43(20) 


Light curve 


6063.5 




17.36(14) 


Light curve 


6064.5 


16.81(01) 


17.41(02) 


Acquisition frame 



Apr 
2009 



Sep 
2011 



CO 



17.0 



17.5 



18.0 



18.5 




5000 



5500 
HJD-2450000 



6000 



Figure 2. Long-term light curve. The circles mark genuine V photometric 
data, while the triangles indicate V magnitudes calculated from measure- 
ments in the R band via the colour index V—R = 0.60 (iPaper ID . The inset 
shows a close-up of the 2012 March and May data. 



an ~1 mag fainter brightness state. This bears consequences for 
the determination of the orbital period. In CVs, a difference in the 
state of bri ghtness can correspond to different radii of the accretion 
disc re.g.. lSmakl[T984l) . potentially affecting the position and the 
(relative) contribution of especially the hot spot. This, in turn, can 
have an effect on the phasing and the amplitudes of radial velocity 
and light curves. 

Table |2] summarizes the photometric data from several observ- 
ing runs. It contains the Heliocentric Julian Day (HJD) correspond- 
ing to the night of the observations, and photometric R and V mag- 
nitudes. The first data point was taken from the calibrated multi- 
colour photometry reported in IPaper l[ which furthermore yielded 
a colour index V — R = 0.60(01) mag. The acquisition frame data 
correspond to the spectroscopic observing runs and were taken in 
the R passband. The respective V magnitude was calculated ac- 

Table 3. Eclipse timings. 



to 


Cycle 


Data 


O-C 


O-C 


(HJD-245 6000) 






(10-4 d) 


(orbits) 


13.8704(23) 


-14 


#20 


2.37 


0.0017 


14.837 88(57) 


-7 


#20 


-6.62 


-0.0048 


15.807 30(23) 





V 


3.80 


0.0027 


19.81824(31) 


29 


V 


-5.34 


-0.0039 


63.672 88(58) 


346 


V 


4.08 


0.0029 


64.6404(21) 


353 


#7 


-4.43 


-0.0032 


64.7775(16) 


354 


#16 


-16.94 


-0.012 



cording to the above colour index. While the colour is supposed to 
change with the brightness of the system we can expect that change 
to amount to much less than 0.5 mag (see e.g. the series of papers 
bv lSpogli. Fiorucci. & Raimond o 2000a. b e). All data points were 
examined for their corresponding orbital phase, using equation ([T]). 
We find that none of them is affected by the eclipse. Specifically, 
the UBVR photometry was taken during phases 0.49-0.77, the ac- 
quisition frame of the low-resolution spectrum at phase 0.66, and 
the other acquisition frames just before the respective spectroscopic 
time series shown in Fig. [T] Lastly, the light-curve data are repre- 
sented as nightly averages of the out-of-eclipse brightness. We can 
thus assume the influence of the orbital variability on the points 
shown in Fig. [2] to be minimal. 

We can clearly see that V728 Sco appears to change frequently 
between brightness states. The UBVR photometry from 2009 May 
20 and the time series photometry from 2012 April 1 find the sys- 
tem at y ^ 18.5 out of eclipse, while the other light curves and 
the time series spectroscopic data catch it at an ^ 1 — 1.5 mag 
higher brightness. Furthermore, the low-resolution spectroscopy 
from 2011 June 29 apparently was taken when the system was in 
between the bright and the low state. While further implications of 
this behaviour are discussed in Section |4l at this point we regard 
the consequences with respect to the determination of the orbital 
period. The light curves in Fig. [T] indicate that the contribution of 
the asymmetrically placed hot spot is strongest in the low state and 
much fainter otherwise. We therefore assume that the eclipsed light 
in high state corresponds mainly to the accretion disc. The latter 
is situated symmetrically around the white dwarf, and thus mid- 
eclipse should represent the superior conjunction of the primary. In 
comparison, in low state the eclipse is much more structured, as 
different sources - the disc, the central object, the hot spot - are 
eclipsed, and the minimum is likely to be offset in phase with re- 
spect to superior conjunction. We thus determined the eclipse tim- 
ings for the high- state data by fitting a second-order polynomial to 
the lowest points of the eclipses, while mid-eclipse for the low-state 
data was calculated from the ingress and egress times of the cen- 
tral object (see Section [33]). These data are summarized in Table 
[3] Here, the fluxes of the spectroscopic data were calculated as the 
total flux in the spectral ranges 6200-7000 A, 4000-5000 A and 
6200-8000 A, for grisms #20, #7 and #16, respectively. They were 
subsequently normalized by dividing through the mean values. We 
then fitted the eclipse timings with a linear regression, weighting 
the data points according to their respective uncertainties. In this 
way we find the ephemeris 



To(HJD) = 245 6015.806 92(19) + 0.138 340(02) E, 



(1) 



where E is the cycle number. 

We would like to point out that no two of our recorded eclipses 
are really the same. Even the two F-band high- state eclipses have 
different depth and show different minor asymmetries and bumps. 
We have tested the robustness of our result by using different sets 
of eclipses for the determination of the orbital period. All permuta- 
tions yield the same result within the errors. 

3.2 Eclipse spectra 

In Fig. [3] we have collected the spectral time-series throughout 
eclipse. These consist of the grism #20 data that show the behaviour 
of the Ha line on two subsequent nights during decline from an out- 
burst, and the two low-resolution series of the blue and the red spec- 
tral range, taken in one night with grisms #7 and #16, respectively. 
The long-term light curve in Fig.|2]suggests that also the latter data 
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Figure 3. Time-series spectroscopic data around the eclipse phases. From left to right: grism #20 data from March 26, from March 27, grism #7 data, and 
grism #16 data, both from May 16. The spectra have been normalised with respect to the continuum and shifted vertically according to their corresponding 
orbital phase. The grism #7 and #16 spectra have been scaled individually to better visualise the sequence. 



were taken during decline, at an only slightly lower brightness state 
than the grism #20 spectra. 

The high-resolution data (grism #20, the two left-hand plots) 
clearly show a double-peaked line profile out of eclipse. Through- 
out eclipse (phases 0.95-1.05) the line presents a single peak, 
which rapidly changes its position in mid-eclipse. This is the well- 
known "Z wave" behaviour, as first the part of the accretion disc 
with velocities towards the observer (that form the blue peak of the 
double-peak profile) and then the part with velocities away from 
the observer (the red peak) are eclipsed. This proves that V728 Sco 
is a disc system. 

The grism #7 and #16 data have too low S/N to analyse the 
mid-eclipse spectra in detail, but they still present a couple of inter- 
esting phenomena. First, it can be clearly seen that the Hell A4686 
line that out of eclipse rivals the H/3 emission line in strength disap- 
pears completely within the limits of the S/N during mid-eclipse, 
while H/3 remains visible (third plot from the left). This presents 
support for the existence of a hot inner disc (as discussed in Sec- 
tion [33]). Secondly, we note that the slope of the red part of the 
spectrum (right plot) in mid-eclipse switches from negative to pos- 
itive towards longer wavelengths. While such behaviour is expected 
during the eclipse of the hotter parts of the system, the contin- 
uum also becomes significantly more "bumpy" indicating that at 
this point and wavelength range the secondary star is the dominant 
light source. For that reason, there is a good chance that a high-S/N 
spectrum taken in mid-eclipse would be able to reveal the spectral 
type of that component. 

3.3 Radial velocities and Doppler tomography 

The Ha emission line in our data presents a broad and double- 
peaked line profile. We thus measure the radial velocities using 
the d ouble-Gaussian technique dSchneider & Youndll98Ql : lshaftej 
Il983h . Here, the line profile is convolved with two identical Gaus- 
sian curves, whose centres are separated along the wavelength axis 
by a quantity d. The wavelength where the fluxes under the two 
curves are identical yields the Doppler shift of the line and thus its 
radial velocity. By varying d and by employing sufficiently nar- 
row Gaussian functions, different parts of the emission line can 
be measured. The usual assumption is that the centre of the emis- 



Table 4. Possible choices for the radial velocity parameters. 



7 




^0 


d 


Criterion 


(kms-i) 


(kms-i) 


(orbits) 


(A) 




-113(6) 


117(7) 


-0.031(11) 


28 


^0 


-92(9) 


128(3) 


-0.037(03) 


32-44 


Ki ~ constant 


-46(7) 


150(8) 


-0.075(11) 


52 


Noise 



sion line represents the part most likely affected by potentially 
present isolated emission sources (e.g., the hot spot or the sec- 
ondary star), while the line wings are mostly formed in the inner 
disc and hence more faithfully reflect the motion of the primary 
star (see iHome & M^ll986h . 

The useful lower limit for the FWHM of the two Gaussians 
is given by the spectral resolution of the data. For the grism #20 
set we therefore use Gaussians with FWHM = 4 A. As mentioned 
above, our two data sets were taken on two subsequent nights dur- 
ing decline from outburst. While the system's brightness during that 
time changed by only slightly more than 0. 1 mag, a visual inspec- 
tion (Fig. [3]) reveals significant differences both in the strength of 
the line and in the line profile. This is corroborated by comparing 
the equivalent widths and radial velocities for corresponding or- 
bital phases. We are thus not allowed to mix the two data sets, and 
since the March 26 data do not cover a complete orbit, we choose 
the March 27 data for the further analysis, excluding the spectra 
around the eclipse (orbital phases 0.9-1.1). 

The radial velocities measured in the way outlined above are 
fitted with a sinusoidal function of the form 

Vr{(p) = 7 - i^i sm{(p - (fo), (2) 

where 7 is the constant term, Ki the semi- amplitude and cp the 
orbital phase with respect to (po. The phases cp have been calcu- 
lated according to equation ([T]), and so one would expect cpo = 
for equation ^ describing the motion of the white dwarf with cpo 
marking the point of its superior conjunction and thus coinciding 
with the photometric eclipse. In order to identi fy the region of the 
line that follows the motion of the white dwarf, Ishaf tej (1 1 9 8 3h in- 
troduced the so-called diagnostic diagram, where the parameters 



Life after eruption. II. 5 



C/) 

E 



E 



100 


-100 

-200 
1501- 

100 

50 





E 



CO 

H — ' 

o 



20 

10 - 


0.4 

0.2 

0.0 - 
-0.2 - 
-0.4 



i h 



— r-p-^ — r-p — ^ 
I I 
I I 
I I 
I I 

I ' . ♦ 

I I 
I I 
i 



1 



H [ 



^ — h 



*. — 



\ 4 
"1 



I 
I 
I 
I 
I 
I 
I 

o 
I 



o o 



\ — h 



\ — h 



• T 

I 
I 
I 

J Ll 







20 



40 



60 



d[A] 



Figure 4. Diagnostic diagram for the March 27 data. The dashed horizontal 
Hnes in the 7 and Ki plots mark the adopted average values for the range of 
separations (i = 32 to 44 A (dashed vertical lines). In the plot of the orbital 
phase (bottom) it instead marks the zero point as determined by Eq. {T}. 



of equation ^ are plotted as a function of the separation d of the 
two Gaussians. The diagram for our data is shown in Fig. |4] The 
parameters present strong fluctuations at small separations corre- 
sponding to the centre of the line, remain comparatively constant 
at medium separations and become variable again at larger separa- 
tions when the noise begins to dominate the line wings. The plot 
gives us three possibilities to choose the parameters, summarized 
in Table |4] First, we can use the point where cpo is closest to 0, cor- 
responding to photometric eclipse. Ignoring the disturbed central 
part of the line this is the case for a separation d = 28 A. Secondly, 
we can take the average for the range d = 32 — 44 A, where Ki is 
constant. And thirdly, d = 52 A represents the last point before the 
error of the semi- amplitude a{Ki) increases sharply, thus indicat- 
ing the start of the noise-dominated parts. Since a closer inspection 
shows that Gaussians with d = 28 A still sample part of the peaks 
of the line, and since the radial velocity curve corresponding to 
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Figure 5. Doppler map of the Ha emission line of the March 27 data 
set. The lower plots present the original (left) and the reconstructed (right) 
trailed spectra. The intensity is scaled to show the range of 1% to 90% of 
the maximum value. 



d = 52 A already presents a large amount of noise, we favour the 
second criterion and adopt the respective parameters for the further 
analysis. 

Fig. \5\ sh ows th e Doppler map for the same data set using 
the code from ISpruii (1998) with th e IDL environment being re- 
placed with a ESO-MIDAS routine dTappert et all l2003h . The in- 
put data consist of the phase- sorted spectra (again excluding the 
phases around eclipse) and the 7 velocity. The latter is iteratively 
adjusted by comparing the average input spectrum with the recon- 
structed one. The final result was 7 = —107 kms~^. Doppler 
mapping visualizes the emission distribution in the system in ve- 
locity space. The orientation of the map is such that it shows the 
system at orbital phase 0.25, so that the primary star has veloc- 
ities (vx^Vy) = (0, —Kwd) and the secondary star accordingly 
(0, i^Ro), where Kwd and K-rd are the semi- amplitudes of the 
radial velocity curves of the white dwarf and the red dwarf, re- 
spectively. The centre-of-mass is situated at (0, 0) and is marked 
in the plot as the point of intersection of the vertical and the hor- 
izontal line. The two lower plots show the trailed spectrogram (in 
phase) from the input data (left) and the one reconstructed from 
the Doppler map (right). We see that the velocity behaviour of 
the emission components is well reproduced; however, the varia- 
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tion in intensity is not. This is due to Doppler tomography assum- 
ing that all emission so urces are constant in intensity at all phases 
(lMarsh&Homel ll988V This formal prerequisite is in practice al- 
most never fulfilled, and especially not in high-inclination systems, 
where the individual emission sources at certain phases are ob- 
scured by other components. However, in the Doppler map this 
should only affect the relative intensities of the emission sources 
but not their positions. 

We then find a prominent isolated component at velocities 
{vx^Vy) = (-25,191) km s~^. Since the resolution of the map is 
50 km s~^ per pixel, the offset in Vx is sufficiently small to iden- 
tify the orbital phase of the emission with that of the secondary 
star. For a more precise measurement we have fitted the component 
with a Gaussian in those spectra in which it was clearly visible. 
A Fourier fit to these data yielded i^add = 221(7) km and 
(fo = 0.52(1) orbits, which confirms above interpretation. Apart 
from the emission lines our grism #16 spectra do not show any ev- 
idence of spectral features from the secondary out of eclipse, and it 
thus can be assumed that the observed emission component is not 
due to magnetic activity, but caused by illumination from the white 
dwarf or the accretion disc. The emitting source will then be lim- 
ited to the side of the secondary star facing the white dwarf, and its 
velocity will not correspond to the centre-of-mass of the secondary 
star, but to a point in between that and the Lagrangian Point Li . We 
will come back to that in Section [3l4l 

The area of the isolated emission is somewhat extended to- 
wards lower Vy and more negative Vx velocities, and is likely to 
represent emission from the gas stream or close to the hot spot. 
The latter itself is optically thick and does not present any (sig- 
nificant) contribution to the Ha line. This can also be seen in the 
behaviour of the Ha flux (lower plot of Fig. [6]) which is declining 
as soon as the hot spot comes into the line of sight. We also see 
that this flux forms a second local minimum at orbital phase 0.45. 
A possible explanation is that part of the Ha emission is produced 
in the gas stream from the secondary and/or a shock front close 
to the hot spot, and that this part is eclipsed by the hot inner disc 
or the optically thick hot spot (see the sketch in the top panel of 
Fig. (6]). Emission from the accretion disc is present, although its 
shape is not that of a closed circle but is rather more reminiscent 
of that of spiral waves (e.g., Steeghs 2001). This would not be an 
unexpected feature in an outbursting CV, but here might just be an 
artifact, since the present data are very limited in both time- and 
spectral resolution, as well as in S/N. 

3.4 System parameters from the spectroscopy 

In order to obtain a more detailed picture of V728 Sco we need 
to know its system parameters. Since the present data lack com- 
pleteness and quality we have to make several assumptions in the 
process. Consequently, this section will not describe the definite 
but a possible configuration of V728 Sco. Our assumptions are the 
following: 

(i) The Ki value extracted from the diagnostic diagram reflects 
the motion of the white dwarf. While the double Gaussians did not 
sample the extreme wings, the constancy of Ki over a wide range 
of separations makes it a good candidate for i^wD- 

(ii) The additional emission component corresponds to an illu- 
minated surface on the secondary star. The middle plot of Fig.[6l 
shows a close-up of the photometric light curves to emphasize the 
behaviour out of eclipse. The signature of the illuminated surface of 
the secondary star can be clearly seen as a broad maximum opposite 
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Figure 6. Top: Sketch of the principal elements of V728 Sco. The orienta- 
tion of the system with respect to the observer is indicated for orbital phases 
0.0, 0.25, 0.5, and 0.75. Middle: The out-of-eclipse light curves for 2012 
March 28 (circles) and April 1 (stars). Both curves have been displaced ver- 
tically to fit into the plot. Bottom: Variation of the flux under the Ha line 
throughout the orbit. 

the eclipse phase. While we do not know the extent of that surface, 
it is safe to assume that its velocity K2 lies in between the velocity 
of the first Lagrange point i^Li and the velocity of the centre-of- 
mass of the red dwarf K-r-d, and it probably will be closer to the 
former of those two. At the same time, the secondary star is elon- 
gated due to the Roche lobe deformation, and in the orbital plane 
the distance from Li to the secondary's centre-of-mass is larger 
than the equivalent radius of the red dwarf. For the sake of simplic- 
ity, we will assume that the distance aRD from the centre-of-mass 
of the secondary to the centre-of-mass of the system is the sum of 
the distance of the centre of the additional emission component to 
the latter and the radius of the secondary, i.e. ord ftadd + Rrd- 

(iii) The system is seen at a high orbital inclination. This is less 
an assumption rather than a conclusion from the presence of a deep 
and structured photometric eclipse. We have performed the subse- 
quent calculations for several inclinations i ^ 75°, finding only 
small differences in the respective results. Taking into account the 
photometric results (Section l33t we adopt i = 82°. 

(iv) Our final assumption regards the type, size and mass of the 
secondary star. It is the least certain of al l, because it is mainly 
based on statistic s. 'Beuermann et al.l (Il998 h find, based on the cat- 
alogue of iRitter & Kolb ( 2003), that the secondaries of the CVs in 
the period range of 3-4 h have spectr al types M3-M5 , with most of 
them being clustered at M4V. Bovaiian et al ](|20l2) conducted an 
interferometric survey on late-type dwarfs. Their result for spectral 
type M4 is based on a single object, GJ 699, for which they deter- 
mine = 0.19 R© and M = 0.15 M©. In CVs, the respective 
parameters are likely to be differen t. A linear-weighted fit to the 
secondary masses collected in the Ritter & Kolbl (l2003h catalogue 
(version 3.18) for CVs with orbital periods between 3 and 4 h yields 
the formal relation 

Mrd = -0.125(8) + 3.03(51)Porb (3) 

for Mrd in Mq and Porb in d. For the fit we have excluded 
masses without error estimates, as well as three 'outliers' (UU 
Aqr, CN Ori, IP Peg). For the orbital period of V728 Sco this 
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Table 5. Parameters from on the eclipse analysis according to lRittej Il980l) 
and the Roche lobe geometry for different orbital inclinations i. 



i 




Gs/a 


Rd/a 


a 


q 


n 








n 




90 


0.256 


0.819 


0.325 


48 


0.21 


85 


0.271 


0.834 


0.318 


50 


0.25 


82 


0.291 


0.852 


0.312 


54 


0.36 


80 


0.309 


0.866 


0.308 


56 


0.45 


75 


0.364 


0.898 


0.300 


62 


0.83 



yields Mrd = 0.29(5) M©. Simply computing the average for 
the three CVs in the immediate period neighbourhood (DW UMa, 
V603 Aql, V1315 Aql, with orbital periods between 3.27 and 3.36 
h) yields Mrd = 0.297(6) M©, proving that above fit is reason- 
able. We then can calcul ate the radius using the equation (16) from 
iKnigge, Baraffe, & Patterson C2011.) that corresponds to our mass 
and period regime, 

Mrd x0.69(3) 



i?RD = 0.293(10) ( 
obtaining i?RD 



-r 



O.2M0 
0.38(14) R0. 



(4) 



With these assumptions, we can calculate the velocity of the white 
dwarf 



'^^WD — Ki I sin i 

and its distance to the centre-of-mass 



ttWD — 



27r 



The corresponding distance for the secondary star is 

. -Porb-^add . 7-, 
ttRD — Cladd + ^RD — — : : V ^RD , 

27r sm % 

yielding the distance between the stellar components to 

a — awD + ttRD • 

The velocity of the red dwarf calculates as 

VRD = 27raRD/Porb , 

and the mass ratio results to 

q — ^;wD/'i^RD, 

which implies a white dwarf mass of 

MwD Mrd/^ . 



(5) 



(6) 



(7) 



(8) 



(9) 



(10) 



(11) 



All these parameters are only weakly dependent on the incli- 
nation i. The parameters (|5} and (|6j can directly be derived from 
the observed quantities, yielding the ranges 128.5 < ^;wd < 
132.5 km s"^ and 0.35 < awD < 0.36 R© for 85° > i > 75°. 
All the remaining parameters depend on the adopted value of Mrd 
and the respective i?RD from Eq. (|4]). We determine in ([Tol a 
mass ratio q — 0.36, and according to fTTI ) a white dwarf mass 
MwD = 0.81 M0, a value very typical for other CVs. We will 
adopt these values for the further analysis of the system parame- 
ters. 



3.5 Analysis of the low-state eclipse 

As shown in Fig. [T] we observed V728 Sco on 2012 April 1 in 
a brightness state that is about 1.5 mag fainter than the brightest 
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Figure 7. The ingress and egress times of the central object and the hot spot 
as seen in the low-state eclipse. The times indicated by the dashed lines 
are not measured, but implied by the width of the corresponding opposite 
transition. 



Table 6. Physical parameter estimates for V728 Sco. 



i 


82° 




0.36 


MwD 


0.81 M0 


Mrd 


0.29 M0 


^RD 


0.38 R0 


a 


1.31 R0 


as 


1.12 R0 


Rd 


0.41 R0 


Rco 


0.09 R0 




0.10 R0 



observed value for the post-nova (from 2012 March 25). There are 
also remarkable differences in the corresponding light curve, show- 
ing a pronounced hump maximum at phase 0.85 and an eclipse light 
curve with a clear distinction of two totally eclipsed, but phase- 
shifted bodies: first that of a luminous source around the white 
dwarf location, which we will call 'central object' (suffix CO), and 
later that of the hot spot (HS) (Fig. Ui- Unfortunately, due to the 
very faint flux at eclipse minimum (V^ ^ 21 mag) some of the 
data points around these phases have rather large errors. In spite of 
this, we here present an attempt to determine some of the binary 
parameters in a preliminary manner. 

For thi s purpose, we used the geometrical method described 
in detail by iRittej (ll98Qh for the analysis of the data on Z Cha, 
a dwarf nova with similar eclipse characteristics. From our light 
curve, we derived for V728 Sco a total mean duration of 297 s for 
the ingress and egress of the central object, and 378 s for those of 
the hot spot. The total eclipse durations (the differences between 
the mean ingress and egress epochs) are 974 s for the central object 
and 1188 s for the hot spot. Finally, the time difference between 
both mid-eclipses was measured to 568 s. Applying equations 8, 9, 
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10 and 11 of lRitteil(ll98Q h we can derive several relations between 
binary parameters for different orbital inclinations, in particular the 
radius of the secondary star Rw, the binary separation a, the dis- 
tance between the secondary star and the hot spot as, the disc radius 
Rd as the distance between the white dwarf and the hot spot, and 
the angle a between the hot spot and the secondary star and with 
the compact object as vertex. If we adopt that Rw is equal to the 
volume radius of the Roche lobe filli ng secondary, we can apply 
equations (2.5b) and (2.5c) in lWarneil (ll995) in order to derive the 
mass ratio q = Mrd/Mwd from i^Ro/a. The results are listed 
in Table |5] Apparently, the mass ratio q is strongly dependent on 
i in the present analysis. On the other hand, there is only a weak 
dependence between q and i in our spectroscopic analysis (Section 
13.41 ). We therefore use the spectroscopically determined mass ratio 
q = 0.36, in order to find the most probable orbital inclination i 
for which this mass ratio is valid, according to the eclipse analysis. 
This yields i = 82°. 

We can now perform an approximate check on the consistency 
of our results. Using i?RD = 0.38 R© from Eq. (|4]), we derive 
a — 1.31 R© from the photometry. This fits well with the spec- 
troscopically determined a — 1.34 R©. One has to keep in mind, 
however, that those results are not completely independent, because 
they both depend on our choice for i^RD- Another test can be done 
by using Mrd, Mwd, and Porb to calculate a from Kepler's Third 
Law, 



G'(MwD + MRD)P,2rb 



47r2 



(12) 



where G is the gravitational constant. This results in a = 1.16 R©. 
This differs from the photometric and spectroscopic results by 1 1 
and 13 per cent, respectively, which still appears acceptable, re- 
garding the uncertainties involved in our calculations of awD and 
aw. 

In Table [6l we present a first parameter estimation of V728 
Sco, using the photometrically determined value for the binary sep- 
aration a. The radii of the central object and the hot spot {Rco 
and i ^HS, respectiv ely) were calculated using equations 16 and 17 
from iRitted (Il980h . The radius Rco is at least a factor 6 larger 
than a typical white dwarf radius, an important observational fact 
which will be discussed in Section ID The observed hump maxi- 
mum at phase 0.85 mentioned above would correspond to an angle 
a — 54° , in agreement with the calculated value in Table [5] This 
can be considered as an additional argument in favour of an orbital 
inclination near 82°. Finally, we can compare the spectroscopically 
determined distance from the white dwarf to the centre-of-mass 
a WD = 0.36 R© with the photometrically derived radius of the 
accretion disc Rd = 0.41 R©. This yields the distance from the hot 
spot to the centre-of-mass to ans ^ 0.05 R©, which corresponds 
to a comparatively small Keplerian velocity ^;hs = 18 kms~^. 
This fits quite well with our interpretation of the elongated exten- 
sion of the emission from the secondary star seen in the —Vx^ +Vy 
quadrant of the Doppler map (Fig.|5}. Note that while the latter was 
taken during decline from outburst, the photometric analysis corre- 
sponds to the low state, a nd thus the disc radii can be expected to 
be different (iSmakll 19841) . Still, due to the apparent short outburst 
recurrence time (Fig. |2) the changes in the disc radius during an 
outburst cycle are probably smaller than in other CVs. 

Since there is only a single eclipse light curve available for our 
analysis, it was not possible to give any error bar on the estimates 
presented here. Furthermore it is worth repeating at this point that 
several parameters depend strongly on the mass of the secondary 
star. Its assumed value is exclusively based on statistics. The pa- 



rameters in Table |6] should therefore not be taken at face value, and 
will have to be tested by future observations. 



4 DISCUSSION 

It is of special interest that V728 Sco shows strong evidence of 
having a comparatively low mass-transfer rate. Its orbital period 
Porb = 3.32 h places it right into the regime of the SW Sex stars. 
These CVs are a g roup of nova-likes that dom inate the orbital pe- 
riod range 2.8-4 h teodrig uez-Gil et al.""2007b). The typical defin- 
ing S W Sex features suggest that these are high mass tran sfer sys- 
tems jRodrfguez-GiL Schmidtobreick. & GansickelE007ah . This is 
supported by their white-dwarf temperatures which is an indica- 
tion for the accretion rate and which exceeds the values expected 
from angular momentum loss through the standard magnetic brak- 
ing scenario (Townsley & Gansicke 2009). Recent investigations 
have shown that most of the non-magnetic CVs in this range and in 
fact all nova-likes se em to be of SW Sex type. T his includes, e.g., 
the old nova RR Pic ( Schmidtobreick et ani2003h . It has been sug- 
gested that SW Sex stars represent a stage in the secular evolution 
of CVs and that, in general, CVs reaching that period range share 
their characteristics (Schmidtobreick et al., in preparation). 

In addition, most old novae ar e indeed high mass-transfer 
systems dlben. Fuiimoto. & MacDon ald 1992), and so V728 Sco 
appears as something of an oddball. However, it is not alone. 
The old nova XX Tau also sho ws the spectral appearanc e 
of a low mass-transfer syst em (ISchmidtobreick et aP l2005l) . 
'Rodri guez-Gil & TorresI (l2005l) suggest an orbital period close to 
3.26 h, i.e. in the immediate vicinity to V728 Sco. Another ex- 
ample is the system V446 Her, which ha s a longer orbital pe- 
riod (4.94 h; iThorstensen & TavloJ 2000l), and shows 'stunted' 
dwarf-nova-type outbursts (Honevcutt. Robertson. & Tumejl 19951 : 



iHonevcutt. Robertson. Turner. & HendenI Il998al) similarto what 
is implied by the long-term light curve of V728 Sco. 
iHonevcutt. Robertson. & Tumej (Il998bh give a few more exam- 
ples. 

This behaviour of some post-novae has been seen as a strong 
indication for the validity of the hibernation model which predicts 
changes o f the mass-transfer rate in the evolution o f the pre- and 
post-nova dShara et all 1 9861 : IPrialnik & Sharalll986h . In that model 
the secondary star suffers an enhanced irradiation by the hot post- 
eruption white dwarf, is strongly driven out of its thermal equilib- 
rium, and as a consequence maintains a much larger volume than 
according to its mass. The gradual cooling of the white dwarf per- 
mits the secondary to relax into its equilibrium state, ultimately 
losing contact with its Roche lobe, thus terminating the transfer 
of mass via Li, and so the post-novae enters into 'hibernation'. 
Hence, post-novae immediately after eruption are expected to drive 
a very high mass-transfer rate tha t will gradually decrease with time 
(iKovetz Prialnik. & Shara|[l988b . Once the mass-transfer rate has 
crossed the limit that allows for disc instabilities, dwarf -nova- type 
outbursts will occur (e.g., Osaki 2005, for a review). The low-state 
behaviour of V728 Sco could be interpreted in this context, es- 
pecially because we are here confronted with a 'rather-old' nova, 
more than one century after eruption. Perhaps V728 Sco is just en- 
tering its transition from a high mass-transfer post-nova stage to- 
wards the onset of hibernation. 

A different explanation for the ob served large spread in mass- 
transfer rates for CVs was given bv iKing et all (Il995h . In their 
approach, the interplay between irradiation from the white dwarf 
and an increase in the Roche lobe radius during the phases of 
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high mass-transfer causes long-term 10^ — 10^ yr) oscilla- 
tions of the secondary's radius, inducing corresponding variations 
of the mass-transfer rate. Here, the fact that most post-novae seem 
to have high accretion rates would simply be a consequence of 
nova eruptions having a higher probability to occur during phases 
of high mass-transfer. This model predicts identical behaviour of 
the pre- and post-novae (beyond the immediate aftereffects of the 
nova eruption). With V446 Her, there is at least one low mass- 
transfer CV known t hat represents a strong support for this picture 

dCollazzi et a02009h. 

Finally, I Schreiber. Gansicke. & Cannizzol (l200d) explain the 
low amplitudes of the outbursts in post-novae with the system 
maintaining a hot ionized inner disc due to irradiation by the 
eruption-heated white dwarf. This limits the area of the disc that 
can contribute to the outburst. As the white dwarf cools down, 
the size of the hot inner disc will gradually shrink, the out- 
burst amplitude will increase, and the outburst frequency will de- 
crease. The very similar amplitudes of the va riation in pre- and 
post-eruption V446 Her (ICollazzi et al.l l2009l) somewhat contra- 
dict that picture. Furthermore, there is as yet no evidence for 
a change in outburst a mplitude or frequency in that post-nova 
(iHonevcutt, Robertson , & Kafka 2011). On the other hand, our 
analysis of the low state eclipse of V728 Sco yields a radius of the 
central object ^0.09 R©, at least a factor ^6 larger than any white 
dwarf radius within the reasonable white dwarf niass range, and 
even more so for our estimated mass of 0.81 M© JProvencal et al.l 
Il998h . This represents strong evidence for the existence of the pre- 
dicted hot inner disc, which is furthermore supported by the pres- 
ence of the strong HeH A4686 emission line in the spectra which 
disappears during eclipse (Section|3]2j, and thus can be assumed to 
originate in that small part of the disc. Also, the low outburst am- 
plitude of 1.5 mag and the high outburst frequency fit well within 
the framework of the Schreiber et al. model. Future long-term mon- 
itoring is desirable to properly investigate the characteristics of the 
outburst behaviour. 



5 SUMMARY 

We have presented photometric and spectroscopic data on the re- 
cently recovered old nova V728 Sco, 150 yr after its eruption, that 
reveal it as an eclipsing system. Its orbital period Porb = 3.32 h 
places it in the midst of a region in the period distribution of CVs 
that is dominated by high mass-transfer systems. Still, the spectro- 
scopic characteristics and the long-term variability that resembles 
dwarf-nova type outburst behaviour point to a comparatively low 
mass-transfer rate. 

The low- state data are of special interest because its light 
curve is characterized by a total eclipse of two different bodies, that 
of a 'central object' surrounding the location of the white dwarf 
and, phase shifted, that of the hot spot. This behaviour is well 
known for SU UMa type dwarf novae wit h Porb < 2 h as Z Cha 
(IWamei]ll974 ), OY Car dVogt et aLlllQsTI) and HT Cas dPattersonl 
but, to our knowledge, it was never observed in old no- 
vae or any other type of cataclysmic variable above the period gap 
(Porb > 3 h). An inspection of published data of about 30 eclipsing 
old nova remnants revealed always V-shaped eclipse light curves, 
similar to those of V728 Sco at high state (Fig.[T]), sometimes with 
short totality phases near eclipse centre. The analysis of that low 
state eclipse provided strong evidence f or the existence of a hot 
inner disc, as predicted for post-novae bv lSchreiberetal.lfcOOOh . 

V728 Sco could soon play the role of a comer stone in our 



understanding of post-novae and cataclysmic variables in general. 
This, however, will only be possible if more and better photomet- 
ric and spectroscopic data become available, especially during the 
interesting low state. Unfortunately, this is all but trivial due to the 
faintness of this target. 
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